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Abstract Electric currents in solar active regions are thought to provide the energy 
released via magnetic reconnection in solar flares. Vertical electric current densities 
J z at the photosphere may be estimated from vector magnetogram data, subject to 
substantial uncertainties. The values provide boundary conditions for nonlinear force- 
free modelling of active region magnetic fields. A method is presented for estimating 
values of J z taking into account uncertainties in vector magnetogram field values, 
and minimizing J1 across the active region. The method is demonstrated using the 
boundary values of the field for a force-free twisted bipole, with the addition of noise 
at randomly chosen locations. 

Keywords: Active Regions, Magnetic Fields; Magnetic Fields, Models; Magnetic Fields,| 
Corona 


1. Introduction 


The energy stored in active region magnetic fields produces large-scale solar activity, 
in particular solar flares. Flares often occur at sites in active regions overlying neutral 
lines where fields are strongly sheared, suggesting a connection between flares and large 
scale electric currents in the solar corona. 

The most accurate determinations of active region magnetic fields and associated 
currents are based on measurements of th e p olarisation state of magnetically sensi¬ 
tive photospheric lines ( del Toro Iniestal . 2003). The data are used to produce vector 


magnetograms, which are maps of the vector magnetic field B = (B x , B y , B z ) across 
regions on the photosphere (where x and y refer to the heliographic west and north 
directions respectively, and z is the radially outwards direction). The process of going 
from polarisation measurements to vector magnetogram values involves atmospheric 
modelling, the resolution of an intrinsic 180-degree ambiguity in the direction of the 
component of the field transverse to the line of sight, rotation of coordinate systems, 
and often also involves rebinning of the data. The resulting field values are substantially 
uncertain. Contributions to the uncertainty include errors in the original polarisation 
measurements, inaccuracy of the atmospheric model assumed in spectro-polarimetric 
inversion, errors in the ambiguity resolution, and loss of information in rebinning. 

The measurement errors depend on the instrument, and two types of instruments 
are in use. Spectro-polarimeters record the full Stokes profile at many wavelengths 
across a magnetically sensitive spectral line, whilst magnetographs, or filter-based in¬ 
struments, record polarisation signals integrated with respect to wavelength over a line 
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( Landi DeglTnnocenti and Landolfi . 20041) . The Hinode Solar Optical Telesco pe (SOT) 
Focal Plane Package includes a Spectro-Polarimeter ( Tsuneta et all 120081) . and the 
Helioseismic and Magnetic Image r on the Solar Dynamic Observatory (SDO) provides 
filtergrams ( Scherrer et al. . 12012 ). Spectro-polarimeter typically have a high spectral 
resolution and provide very accurate measurements of the polarisation state, but they 
take longer to produce the data for a vector magnetogram and the magnetic field may 
change during the observation interval. Filter instruments produce data for a region on 
the Sun more rapidly, but they lack the spectral resolution and intrinsic high accuracy 
provided by a spectro-polarimeter. 

Given a set of vector magnetic field values, estimates of the vertical (locally radial) 
electric current density at the level of the photosphere (z = 0) may be made using 
Ampere’s law: 


HoJz(x,y,0) = v- 


where we introduce the notation 

dB x 


dy 


A dB V 

and v = -- 


z=0 


dx 


(1) 


( 2 ) 


i=0 


If the uncertainties in the field values are a x and <j y , and centered differencing is used 
to estimate the derivatives, then the uncertainties in the derivatives are a u = a x /\/2h 
and a v = a y /V^h, where h is the grid spacing in the magnetogram in x and y. If the 
spatial scale for variation of the horizontal field is L , then we have u ~ B x /L and 
and By denote characteristic values of the field components. In 


v ~ By/L , where B x 
that case 


L (j x 

hB x 


and 


hB x 


(3) 


Equations © imply that if there are significant fractional uncertainties in the field 
values, then the derivative estimates may be substantially uncertain, in particular for 
high-resolution data (with L> h). The current densities obtained using Equation dTj) 
will be correspondingly uncertain. 

Errors in t he determination of current den s ities fr om v ector magnetogram data wer e 
discussed bv IPevtsov. Canfield, and Metcall (1994) and Leka and Skumanichl ( 19991) . 
who assigned constant uncertainties for field components in the directions along, and 
transverse to, the line of sight. Uncertainty values may be assigned to individual 
vector magnetogram field values based on the quality of the model fit in the spectro- 
polarimetric inversion, and this in formati on i s now routinely supplied with vector mag¬ 
netogram held values (e.g. iHoeksema et all I2014I) . Other s ources o f un certa inty have 

also been discussed, including errors in ambiguity resolution ( Leka et al I 200i , and er- _ 

rors introduced by the rebinning of data in the construction of magnetograms (Leka and Barnei .[ 
2012 !). 

The nonlinear force-free model is oft en used to model the coronal m agnetic held 
based on vector magnetogram data (e.g. IWiegelmann and SakuraiL 20121) . The model 
is dehned by cur l B = aB and divB = 0, where a is the force-free parameter. The 
Grad and R.ubinl ( 1958 ) method of solution of these equations takes as boundary con¬ 
ditions the values of B z at the photosphere, and the values of J z (or the values of a) 
over one polarity of B 7 . Other methods, in particular the optimization meth od (e.g. 
Wheatland. Sturrock. and Roumeliotid . 2000L Wiegelmann and Inhester , 2010l) , use all 
three components of B at the photosphere as boundary conditions. In practice all 
methods modify the boundary data substantially from the values in th e vector magne- 
togram in order to achieve a solution to the force-free model (see e.g. iDe Rosa et all 
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2009[ Wheatla nd and Leka. 2011. Wiegelmann et aL . 2012 ). For the active region mod¬ 
elled by Wheatland and Lekal (2011 ) using the ‘self-consistent’ Grad-Rubin procedure 
( Wheatland and Regnierl 2009 ). the current densities in the magnetogram were altered 
to the extent that the average absolute change in the horizontal field was 170 gauss, 
and the average ratio of the change in the horizontal field to the assigned uncertainty 
in the horizontal field, on a pointwise basis across the magnetogram, was « 9. The 
changes greatly exceeded the nominal uncertainties, which is typical. The ‘prepro¬ 
cessing’ procedure often used with the optimizati on method of nonlinear force-free 
solution (IWiegelmann, Inhester. and Sakurail 120061) alters the transve r se com ponents 
of the held by < 500 gauss in strong held regions ( Fuhrmann et all 12011 ). If non¬ 
linear force-free codes are applied to vector magnetogram boundary data without 
these changes, accurate solutions to the nonlinear force-free model are (generally) not 
obtained. 

The large changes required in vector magnetogram data to achieve nonlinear force- 
free model solu tions are most l i kely due to the inconsistenc y of boundary data with the 
model (see e.g. Schriiver et al . 2008 . De Rosa et ali . 20091) . The solar magnetic held is 
not expected to be force-free at the level of the photosphere. Errors in the held estimates 
also lead to inconsistency. Depending on the solution method, this problem may cause 
the soluti on magnetic helds to have signihcant departures from the divergence-free 
condition ( Valori et al. . 120131) . 

In this paper we reconsider the problem of estimating electric current densities from 
vector magnetogram boundary data, taking into account uncertainties in the data. A 
method is introduced for calculating a set of boundary values of J z which minimize 
departures from centred difference estimates for the current densities at locations where 
the values are accurate (according to the uncertainties), and which minimize the sum 
of J? across the magnetogram. The goal is an estimate of the current density which 
avoids large values produced by errors. The method is tested in application to a test 
case with known errors. The paper is divided as follows. Section 2 presents the method, 
and Section 3 describes the test and analyzes the results. Section 4 draws conclusions. 


2. Method 


Equation © provides the vertical current density in terms of the horizontal field gradi¬ 
ents u and v. The problem consists in estimating the gradients from the data, subject to 
unertainties, and we also seek to avoid large (spurious) values of the resulting vertical 
current densities, produced by points with large uncertainties. Hence we consider the 
problem of minimizing 


F = E 


(Ujj - O 2 
2(j2 uij 


(Vjj - <f ) 2 

2(7 vij 


+ A El Jz ij > 
ij 


( 4 ) 


where the indices refer to points ( Xi , yj) in the magnetogram, u and are estimates 
of the gradients, and where A is a constant. It is useful to non-dimensionalize by 
expressing lengths, field strengths, and current densities in units of characteristic values 
B s , L s , and B s /yoL s , respectively. Using Equation ©, Equation © may be written 
in non-dimensional form 


F = E 


{Uij - O 2 
2a liJ 


(■Vij ~ ffi ) 2 

2<J vij 


+ A E(% 



( 5 ) 
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where A = XBg/fx 0 L^. 

A minimum value of F is obtained when dF/duij = dF/dvij = 0, which gives two 
coupled linear equations for the values of and v-ij corresponding to the minimum: 


..min _ .est 

a ij U ij 


u%j 

..min _ ..est 

u ij ij 


-2A =0 

+ 2A « n -u“ in ) =0. 


( 6 ) 


The simultaneous solution of Equations © is given by 


M min = _ 
« d 


‘ + »A PS 


VIJ 


uij 


?; est 

2A— 


VIJ 


(7) 


and 


, est 

2A—Zr— 


'UIJ 


'UIJ 


\ u est 

2A I —B— 


'VIJ 


( 8 ) 


where 


d = 


2 2 
ij ij 


2A 


2 2 
®u ij ®v ij 


(9) 


The current density corresponding to the gradients which minimize F is J™ n = i;™" — 
which evaluates to 


rest 


Tirnn _ 

zij “ i • 2.w;v 


( 10 ) 


where 


rest _ est _ est 
17 2 ij u ij 11 ij 


(ii) 


and 


2 _ 2 I 2 

ij ® u ij ij ’ 


( 12 ) 


Equation m provides a surprisingly simple solution to the problem: the current 
density at each point is reduced from the usual estimate by a factor 


f ij = (l + 2Aa a JiJ ) \ 


(13) 


so at points with larger uncertainties the current density is reduced more. We refer to 
Equation ca as the ‘minimum current’ estimate for the current density. 

In the following we use centered differences to provide the estimates of the gradients: 


?/ est = 
u ij 


B x (Xj , yj + h) - B x (xj , yj - h) 
2 h 


(14) 


and 


u est = 

u ij 


Byi^Xi -f- h, yj) By (xj h, yj) 


2 h 


(15) 
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where h is the grid spacing in the magnetogram in x and y. With these choices the 
uncertainties in the gradients, using propagation of errors, are 


u ij 


1 

2 h 




2 

x ij +1 


+ C7 


2 

yij —1 



(16) 


and 


ij 


1 

2 h 




2 

y 


+ tT 


2 

yi—ij 


) 


1/2 


and the error in the current density is 


(17) 






2 

x ij -\-1 


+ <7 


2 

yij -1 


+ CT 


2 

2/ *+1j 


+ cr 


2 

yi—lj 



(18) 


Higher order differencing schemes are sometimes used. In the general case Equation m 
is replaced by 

^ H c k B x (xi, yj + kh), (19) 

k 

where the sum enumerates the points involved in the differencing, and the c k are 
coefficients. Equation m is replaced by 


U?/ ij 


1 

2 h 


{ a l a tij+k ) 1/2 . 

k 


( 20 ) 


and corresponding relations hold for vff and a v ij. Equation (fTSl) is replaced by 


° Jij ~ 2 h 


^ ^ OL k {&xij-\-k d" &y ij+k) 
. k 


1/2 


( 21 ) 


3. Test using simulated data 

3.1. Boundary conditions and minimum current estimate 

The method is demonstrated in application to a test case consisting of boundary 
conditions for a twisted bipole, with noise added to the boundary values at some 
points. 

Panels (a), (b), and (c) of Figure 1 illustrate the twisted-bipole boundary conditions. 
Panel (a) shows the vertical component of the boundary field, B z (x,y,0), which is 
constructed by adding the field from monopole sources at locations (x+,y+, 2 :+) and 
(x-, j/_, Z-). In non-dimensional units we have 


B z {x,y,Q) = c[^±- Z -^-y ( 22 ) 

with 

R± = [{x-x±) 2 + (y-y±) 2 + z±] 1/2 , (23) 

and with C chosen such that max \B x {x, y, 0)| = 1. The boundary field is constructed 
over the region 0 < x < 1 and 0 < y < 1, with the choices x+ = 0.6, X- = 0.4, 
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y+ = y- = 0.5, and z+ = Z- = —0.2. A current is introduced by assuming a boundary 
distribution of the force-free parameter a over the positive polarity of the field. We 
choose 


a(x,y,0) 


a 0 if B z > Bi 
0 otherwise, 


(24) 


with Bi = 0.9 and ao = 15. The boundary conditions on B z and a are used to calculate 
a nonlinear force-free field in the cubical regio n defined by the b oundary region and 
0 < z < 1, using the CFIT Grad-Rubin code ( Wheatland . 20071). The calculation is 
performed on a 100 x 100 x 100 grid. Panel (b) of Figure 1 shows the magnitude of the 
horizontal field 

B h (x,y, 0)= [B x (x, y, 0) 2 + B y (x, y, 0) 2 ] 1/2 (25) 


for the calculated force-free field and panel (c) shows the corresponding values of the 
current density J z (x 1 y 1 0) for the force-free field, estimated using Equation ([!]) with 
the centered difference values Equations ED and ED used for the field gradients. 
The patch of positive current matches the current density J z = aB z defined by the 
boundary conditions on B z and a assumed for the force-free calculation [Equations (1221) 
and (1231) ]. and the patch of negative current density follows from the mapping of the 
boundary values of a along field lines in the force-free field. 

Noise is added to the boundary values of B from the force-free calculation, to provide 
a set of ‘observed’ boundary values. At a fraction 6 of the boundary positions, randomly 
chosen, additional components 5B X and 5B y are added to B x (x : y,0) and B y (x,y 1 0) 
respectively, where 5B X and SB y are normally distributed random numbers with mean 
zero and standard deviation a x = cr y = ao, and where ao is a chosen constant. The 
boundary values of B z (x,y,0) are not altered. Panel (d) of Figure 1 illustrates the 
resulting noisy values of the horizontal field, for the choices 6 = 0.1 and <to = 0.25, and 
panel (e) shows the corresponding current density, which we use for J| st . The values 
are obtained using Equation CD applied to the observed (noisy) boundary values with 
centered differencing used to estimate the gradients. The noise in the horizontal field 
leads to spurious values of the current density, which mask the positive and negative 
patches seen in panel (c). 

Panel (f) of Figure 1 shows the minimum current estimate </“ ln for the current 
density calculated using Equation m with the choice A = 5. This value of A was 
chosen because it reduces the noise in weak field regions to small values. The values of 
J z (x,y, 0) shown in panel (e) provide the estimate J| st , and the uncertainties in the 
current density values aj follow from Equations (fl6l) - (H7l) assuming a x = a y = ao at 
points with noise, and a x = a y = 0 at points without noise. Panel (f) shows that J™ ln 
returns a good approximation to the underlying distribution of current density in the 
boundary, except at locations affected by the noise in the horizontal field. 

With the chosen value of A, the RMS value of (Jf st — T“ ln )/crj over the boundary 
points with aj > 0 is close to unity. This is because the current density at the locations 
affected by noise is reduced by a value essentially equal to the noise. For this test case 
the minimum current estimate changes the boundary values within the uncertainties. 


3.2. NLFFF reconstructions 


The minimum current estimates for the current density enable reconstruction of a 
nonlinear force-free field. If the CFIT code is applied using the noisy current density 
estimates J? st the Grad-Rubin iteration does not converge: the noise in the data pre¬ 
vents a successful reconstruction. However, accurate force-free solutions are obtained 
using the minimum current estimate for the current densities. 
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Figure 2 illustrates the force-free solutions. Panel (a) shows the original NLFFF 
obtained from the boundary data without noise. The greyscale shows the boundary 
values of B z , which are shown also in panels (b) and (d), and the white curves are a 
set of field lines originating at the positive pole. (Each of the panels in Figure 2 shows 
a cropped central region, by comparison with Figure 1.) Panel (b) shows the P and N 
solutions obtained from the boundary values of B z together with the minimum current 
method estimates for current density. The white field lines correspond to the P solution, 
and the black held lines to the N solution. Panel (b) shows that the P and N solutions 
are twisted bipole configurations which are similar to the original held [panel (a)] and 
to each other, but are not identical. The P and the N solutions have less twist than the 
original held due to the reductions in current density at points affected by noise, and 
the energies of the helds are correspondingly decreased. The energy E of the original 
held shown in panel (a) is given by E/Eq = 1.083, where Eq is the energy of a potential 
held with the same lower boundary v alues for B z [calculated by the method described 
for the CFIT code ( Wheatlandl . l2007h ]. The energies of the P and N solutions shown in 
panel (b) are given by Ep/Eq = 1.064 and E-^/Eq = 1.048. The reconstructed P and 
N helds contain about 80% and about 60% of the free magnetic energy of the original 
held, respectively. 

We apply also the self-consistency procedure presented in Wheatland and Rearmed 
( 2009h . which uses boundary values at both polarities to identify a single NLFFF 
solution. The procedure involves a cycle of calculations of P and N solutions. After 
a pair of solutions is calculated, the boundary values of a in the two solutions are 
averaged, taking into account uncertainties in the a values. The resulting averaged a 
values provide the boundary conditions for the P and N solutions in the next cycle. 
After a number of cycles, the P and N solutions agree, and provide a ‘self consistent’ 
solution. Panels (c) and (d) in Figure 2 illustrate the results. Panel (c) shows the P 
and N solutions - the white and black curves, respectively - after nine self-consistency 
cycles. The two sets of held lines are mostly overlapping. The energies of the P and N 
solutions are Ep/Eq = E^/Eq = 1.085. A self-consistent solution is found, with a free 
energy a few percent higher than the energy of the original held. Panel (d) shows the 
J z values in the self-consistent held, for comparison with panel (c) of Figure 1. The 
self-consistency procedure recovers a good approximation to the original set of current 
density values, although there are detailed differences. The method is able, to some 
extent, to ‘remove the noise’ because at a boundary point affected by noise, say at the 
P polarity, the uncertainty in the a value is large. The averaging procedure replaces 
this a value with the a value from the N solution, provided it has a much smaller 
uncertainty. This procedure generally recovers the original a values. 

The self-consistent solution is qualitatively simil ar to the origina l held . A quantita¬ 


tive measure is provided by the mean vector error ( Schriiver et al . , 200(1 1: 


Em ~ M ^ 


|B 0 rig ,2 Bsc ,i | 


(26) 


where B or i g is the original NLFFF, B sc is the held obtained by the self-consistency 
procedure, and M = 100 3 is the number of grid points. We hnd E m = 1.6 x 10 -2 , 
indicating close correspondence. 

4. Discussion and conclusions 


A method is presented (the ‘minimum current’ method) for estimating vertical electric 
current densities J z at the photosphere from vector magnetogram data. The procedure 
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Figure 1 . Application of the method to a simple bipole test case: boundary values. Panel (a) shows 
the boundary conditions on B z for the test case and panel (b) shows the horizontal field in the 
boundary, obtained from a NLFFF solution with the chosen boundary conditions on B z and current 
density. Panel (c) shows values of J z in the boundary for the force-free field. Noise is introduced to 
the horizontal field values at a fraction of boundary locations, as shown in panel (d), which influences 
the estimated values for J z [panel (e)]. Panel (f) is the minimum current estimate J| st for the vertical 
current density from the noisy boundary data. 


alters the horizontal field gradients used to calculate the current densities subject to 
the uncertainties in the estimates for the field gradients, and also minimizes the sum 
of the square of the vertical current densities over all points in the magnetogram. 
The resulting minimum current estimate for the vertical current density at a given 
boundary location is reduced from the original estimate by a factor 1 + Aerj, where 
aj is the uncertainty in the original estimate of the current density, and A is a chosen 
constant which determines the relative importance of adherence to the original current 
density estimate or minimization of the sum of J z over all boundary locations. 

The minimum current method provides an estimate for the current density which 
preserves values of the current density at locations with small uncertainties, and re¬ 
duces the current density at locations with large uncertainties. The nonlinear force-free 
(NLFFF) model is often applied to vector magnetogram data to provide a proxy coronal 
field (e.g. IWiegelmann and Sakurai . 20121) . The construction of an accurate solution to 
the model requires modification of the boundary da ta from observed magnetogram 
values (e.g. Schriiver et al. . 20081 De Rosa et al. . 2009). The minimum current method 
may provide a more systematic approach to this modification. 

A demonstration of the method is presented, in application to boundary field values 
for a NLFFF twisted bipole field configuration, with noise of a known amplitude added 
at a fraction of boundary locations. The model field values are taken to represent 
observed data with known uncertainties. The minimum current method recovers a good 
approximation to the original current density values in the boundary from the noisy 
data, except at locations affected by the noise. The method alters the boundary values 
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Figure 2. Application of the method to a simple bipole test case: nonlinear force-free fields. Panel 
(a) shows a set of field lines for the force-free field calculated from the original bipole boundary data, 
originating at the positive pole. Panel (b) shows comparable sets of field lines for the P solution (white 
lines) and for the N solution (black lines), calculated from the bipole values of B z together with the 
minimum current method estimates J™ m for the current density over the P and N polarities of the 
field in the boundary, respectively. Panel (c) shows the solutions obtained with the Wheatland and 
Regnier (2009) self-consistency procedure (the white/black field lines show the P/N solutions at the 
end of a set of self-consistency cycles). Panel (d) shows the boundary values for current density in the 
self-consistency solution. The images in each panel are cropped to a central region in the computational 
domain, to more clearly show the field lines. 


for current density within the uncertainties. The minimum current method bound¬ 
ary values are shown also to enable NLFFF reconstructions. The P and N solutions 
obtained from the boundary values are twisted bipole field configurations similar to 
the original NLFFF, but with reduced energies, due to the reduction in the bound- 
ary current density at location s affected by the noise. The self-consistency procedure 
( Wheatland and Rcgnicil 120091) is applied, and demonstrated to produce a field which is 
a close match to the original NLFFF. The energy of the self consistent solution is within 
a few percent of the energy of the original NLFFF. The self-consistency procedure uses 
information on current density over both polarities of the boundary field, and hence is 
able to recover current density values affected by noise at one or other polarity. 

We have chosen to minimize J z across the vector magnetogram to produce boundary 
data suitable for NLFFF modeling. The vertical current density is a boundary condition 
for the force-free model, and large estimates of J z from vector magnetogram data 
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(which may be due to observational errors) can prevent NLFFF codes from producing 
an accurate force-free solution. However, this choice may remove or reduce local current 
concentrations which are real, and which contribute significantly to the energy of the 
field. Methods of resolving the 180-degree ambiguity us ed to pr oduce vector rnagne- 
tograms may also excl ude real currents in the data (e.g. iMetcall Il994 : Metcalf et al. 


2006 : Leka et all . l2009l) . An advantage of the present method is that, if the uncertainties 


associated with currents are small, the currents are preserved. The accuracy of the 
minimum current estimate depends on the accuracy of the uncertainties provided. 

The test of the method presented here is highly idealised. It is not intended to 
provide a realistic example of the application of the method to vector magnetogram 
data, but only to demonstrate the method, and to show its performance in the best 
case. 

The method of inclusion of noise in the boundary data (adding Gaussian noise of 
known amplitude at specific locations to the horizontal field components of a NLFFF 
calculated from analytic boundary conditions) does not produce a close approxima¬ 
tion to real solar data. We have not (yet) confronted the challenge of estimating 
electric currents from solar data with the minimum current method. Other studies 


have examined the influence o f noise in vector magnetogram data (IWiegelmann et al 


201(1 iLeka and Barnes . 2012 1. These studies add noise to synthetic line profiles, and 


then apply the various steps involved in constructing vector magnetograms (inversion, 
ambiguity resolution, co-ordinate transformations and rebinning, etc.). This process 
accurately mimics the influence of measurement errors on vector magnetograms. Tests 
of the method on data of this kind may be tried in future. 

The current density estimates presented here use centered differences, but higher 
order differencing schemes could in principle be used, following Equations ED (ED. 
High order schemes offer advantage in estimating gradients when applied to well- 
sampled data with small uncertainties, but the advantage is lost in application to 
sparse data, or data with substantial errors, and these schemes have the disadvantage 
of introducing correlations in the neighbouring estimates. As such we have chosen to 
use centered differences. 

In future work the method will also be applied to vector magnetogram boundary 
data. This will allow investigation of whether the method is of use for coronal magnetic 
field modeling. 
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